Back-to-Back Black Holes decay Signature at Neutrino Observatories 
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We propose a decay signature for non-thermal small black holes with masses in the TeV range 
which can be discovered by neutrino observatories. The black holes would result due to the impact 
between ultra high energy neutrinos with nuclei in water or ice and decay instantaneously. They 
could be produced if the Planck scale is in the few TeV region and the highly energetic fluxes are 
large enough. Having masses close to the Planck scale, the typical decay mode for these black 
holes is into two particles emitted back-to-back. For a certain range of angles between the emitted 
particles and the center of mass direction of motion, it is possible for the detectors to measure 
separate muons having specific energies and their trajectories oriented at a large enough angle to 
prove that they are the result of a back-to-back decay event. 



I. INTRODUCTION 

In brane world models with a large extra-dimensional 
volume ^-[3] or in even in four dimensions if there is a 
large hidden sector of particles [4] , quantum gravitational 
effects could become important anywhere between the 
traditional Planck scale, i.e. some 10 16 TeV and a few 
TeV. If the Planck scale, i.e. the energy scale at which 
quantum gravitational effects become important, is in the 
lower end of this energy range, the collision of particles 
can result in the creation of small black holes with TeV 
masses when particles collide with center of mass energies 
larger the Planck scale. 

The formation of black holes in the collision of particles 
has been studied since the 70's. In 1972, K. Thorne pro- 
posed the Hoop conjecture [5] which states that a black 
hole forms whenever the impact parameter b of two col- 
liding objects (of negligible spatial extension) is shorter 
than the radius of the would-be-horizon (roughly, the 
Schwarzschild radius, if angular momentum can be ne- 
glected) corresponding to the total energy M of the sys- 
tem [6j 



While the hoop conjecture is intuitively very satisfac- 
tory, it is not enough to prove that black holes do indeed 
form in such collisions. However, there are now a proofs 
available of the formation of a closed trapped surface 
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when two particles collider at very high energy above the 
Planck scale. The formation of a closed trapped surface 
suffices to demonstrate gravitational collapse and hence 
the formation of a black hole [7] . The proofs [SHTTj cov- 
ers both both zero and non-zero impact parameters. Re- 
markably the proof of Eardley and Giddings is analytical 
in the case of a four dimensional space-time [llj . This 
demonstrates the formation of a classical black hole in 
the collisions with a non-zero impact parameter of two 
particles with energies much larger than the Planck mass. 
This work has been extended to the semi-classical regime 
by Hsu [T2]. Semi-classical black holes are expected to 
have masses in the range from 5 to 20 times the Planck 
scale [15] . 

Most of the articles related to the production of small 
black holes via particle collisions at colliders or in cosmic 
rays have considered semi-classical black holes P3H22] . It 
is however possible that the center of mass energy avail- 
able in such high energy collisions is not large enough to 
create semi-classical black holes. It was therefore pro- 
posed [2"5H2"5] to consider quantum black hole which are 
non-thermal objects with masses close to the Planck mass 
which should be easier to produce. As they are non- 
thermal objects, quantum black holes are expected to 
decay into a small number of particles, typically two. Ex- 
perimental signatures for such decays are very different 
from the one of semi-classical objects which are expected 
to decay into several particles in a final explosion, see e.g. 
[2"61 |2"T] for recent reviews. 

Bounds on the Planck scale using Earth skimming neu- 
trinos creating black holes in the Earth crust have been 
derived in []j|J [201 12H] ■ The back-to-back decay signature 
for quantum black holes produced in cosmic ray events 
was first proposed in [29] . The authors study the possibil- 
ity for the two particle showers produced by the particles 
resulting from the back-to-back decay of the black holes 
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to be spatially separated and detected as two simultane- 
ous shower events by cosmic ray observatories (current 
earth based or future space based experiments). In that 
case study refers to quantum black holes which are gen- 
crated due to the interaction of ultra high energy cosmic 
rays (UHECR) or neutrinos with particles in the upper 
atmosphere and decay instantaneously into two particles 
which move back-to-back in the center of mass reference 
frame. It is shown in that even if a small percentage 
of this type of events can be detected, there is parameter 
space for which detection is possible. 

Here we propose to take the idea one step further by 
analyzing whether it is possible to discover such a black 
hole decay signature in ice or water with the help of neu- 
trino observatories. Besides being produced at colliders 
or in the atmosphere by high energetic collisions of cos- 
mic rays with nuclei, quantum black holes can also be 
produced due to the collision between highly energetic 
neutrinos with nuclei in water or ice. Neutrino observato- 
ries are designed to detect muons induced by high-energy 
neutrinos. Because of the characteristics of the propaga- 
tion of muons in water and ice, the neutrino direction 
can be derived with high accuracy [301 [3T] . As neutri- 
nos only interact weakly and their trajectory points back 
to their sources, energetic neutrino events would point 
directly towards sources capable of producing these en- 
ergetic events. We propose that these observatories can 
be used to search for quantum black hole events. Indeed, 
when neutrinos of high enough energies collide with par- 
ticles in water or ice, if the center of mass energy is larger 
than the Planck mass, quantum black holes can be cre- 
ated. As stated before, black holes with masses close 
to the Planck mass decay preferentially into two parti- 
cles which then produce secondary showers which can be 
seen by the neutrino experiments. 

If the Planck mass is of the order of a few TeV, only 
black holes with masses above this energy scale can be 
produced. This implies that, to form a black hole, the 
energy of the neutrino has to be of the order of 10 7 GeV 
or above. The range of interaction lengths for neutrino 
energies (E v ) between 10 7 - 10 9 GeV is 6.6 x 10 3 9.4 x 10 2 
km water equivalent in rock |32j . This means that the 
Earth is opaque for neutrinos with energies in this range 
or larger. Only Earth skimming neutrinos 33: and those 
coming from above the horizon are thus useful for our 
considerations. 
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II. BLACK HOLES PRODUCTION 



In this section we briefly describe the production cross 
section for quantum black hole formation. Note that the 
formulas are extrapolated from the semi-classical regime. 
The black hole production cross section as a result of a 
neutrino interacting with nucleon [y N — > BH) is given 



a(E v ,x min ,M D ) = J 2zdz 



dxF(n) (2) 



,M D )Y,Mx,Q). 



In this equation Mo is the 4 + n dimensional reduced 
Planck mass, z — b/b max with b the impact parameter 
and b ma x the maximum value of the impact parameter for 
which black hole creation can occur as a result of the col- 
lision between the two particles, x m i n — MBH,min/Mi) 
and n is the number of extra-dimensions. F(n) and y(z) 
are the factors introduced by Eardley and Giddings jTT] 
and by Yoshino and Nambu [3"T| . The Schwarzschild ra- 
dius in 4 + n dimensions is given by 



r s (us,n,M D ) 



us/M, 



D 



1/(1- 



where 



k(n) 



.„- 3 r((3 + n)/2)l 1/(1+n) 



(3) 



(4) 



Furthermore, note that s = 2xm^E Ul with tun the nu- 
clei mass and E v the neutrino energy. The functions 
fi(x, Q) are the parton distribution functions. Black hole 
production by cosmic neutrinos might be suppressed in 
comparison to the production rate from UHECRs |35j . 
However this is a model dependent question. 

Only the flux of highly energetic neutrinos is relevant 
for the present case of study, flux which can be esti- 
mated by considering two sites of productions: at the 
source and between the source and the detection place, 
usually Earth. The production sites of the extragalactic 
UHECR, which include (AGN) and (GRBs), are also as- 
sociated with the ones for neutrinos which are produced 
through pion decay in proton-proton or proton-photon 
interactions within the source [36] and the flux depends 
on the composition of the cosmic rays at high energies, 
which can be pure protons, neutrons, heavy nuclei or a 
mix of these [371138]. 

In order to estimate the number of quantum black hole 
events expected at a neutrino detector like IceCube, we 
use two models for the energy flux of the neutrino at high 
energy proposed in [39]. Here they use a a smoothly- 
broken power law to estimate the fluxes of high energy 
neutrinos reaching the Earth that also include the recent 
observation of two PeV-energy shower events by IceCube. 
These two models are based on a E~ 2 accelerated pro- 
ton spectrum, |40j , and use first a 7r + only decay channel 
and then tt ± and ^ decay channels to produce neutri- 
nos. We combine this with the geometrical acceptance 
of the IceCUBE detector, [41] to find the numbers of the 
quantum black holes produced by a flux of neutrinos in 
the energy range between 10 7 to 10 8 GeV at IceCUBE 
in one year. Depending on the model (i.e. the number of 
space-time dimensions, (0, 1, 2, 3, 4, 5, 6, 7) where corre- 
sponds to the large hidden sector model, n = 1 to Ran- 
dall Sundrum brane world model and higher n models 
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to ADD brane world models), we find (2092, 56342, 2.2 x 
10 5 ,4.9 x 10 5 ,8.5 x 10 5 ,1.2 x 10 6 ,1.7 x 10 6 ,2.2 x 10 6 ) 
black holes per year for the first model of neutrino flux 
and (262,5843,21872,47260,79754, 1.1 x 10 5 , 1.5 x 10 5 , 
2.0 x 10 5 ) for the second one. Note that we have set 
F(n) = 1 and y(z) = 1 to produce our numbers. These 
factors are not known for non-thermal quantum black 
holes but they are expected to be order unity. 

Following [42 , it is easy to estimate the branching ra- 
tios for the decomposition of non-thermal black holes. 
Assuming flavor conservation and baryon number conser- 
vation we find that black holes produced from a u-quark 
and a muon- neutrino will decay 33.3% of the time back 
to + u, [i + d and gluon. The same applies to 
black holes formed by the collision of a e?-quark and a 
muon-neutrino. 



III. BACK-TO-BACK BLACK HOLES DECAY 
SIGNATURE IN WATER AND ICE 

A black hole produced as a result of the collision be- 
tween a highly energetic neutrino of energy E v (we shall 
neglect neutrino masses) and a particle of mass m has a 
rest mass Mbh equal to 



M. 



BH 



Vm 2 +2E U 



m 



and is moving relativistically with 



Ibh 



M 



BH 



(5) 



(6) 



Obviously an on-shell black hole can form only when 
Mbh > Mpi. In the remaining of this paper we shall 
assume that this condition is satisfied. Furthermore, we 
use Mbh and jbh to refer to the black hole mass and 
to the Lorentz factor. Note that the quantum black hole 
mass is determined by the neutrino energy and is a con- 
tinuous quantity. For our calculation and our numerical 
simulations knowing the orders of magnitude for Mbh 
and jbh is sufficient. In principle one could take take 
into account the amount of energy which is radiated via 
gravitational radiation, the dependence of the horizon 
formation and black hole mass on the value of the im- 
pact parameter b and so on. However, assuming that 
the impact parameter is small enough for a black hole 
to form, the mass and Lorentz factor of the black hole 
which forms can vary by less than an order of magnitude 
when considering all the above mentioned effects. 

Quantum black holes are non-thermal and they are 
expected to decay to a small number of particles, most 
likely two. In the center of mass the two particles pro- 
duced by the instantaneous decay of such a black hole 
are emitted back-to-back due to momentum conserva- 
tion. There is no preferred direction in which the decay 
takes place, since the differential cross section is angular- 
independent and the energy of the resulting decay prod- 
ucts is restricted by energy conservation which translates 
into fact that the sum of the two masses m a and mi, has 



to be smaller than Mbh- In the center of mass reference 
frame the momenta of the two particles are opposite vec- 
tors with magnitudes equal to 



[(M 



P= 2M^ - (7) 

The energies and momenta of the two particles in the 
laboratory reference frame (Earth reference frame) are 



(m - m b ) 



or 



?4 



Ibh 
-PbhIbh 



-PbhIbh 
Ibh 



E, 



(8) 



Pi± = Pi-L 



where i = a, 6; Ei and pt are the energy and momen- 
tum for the i-th particle measured in the center of mass, 
while the primed quantities are the corresponding ones 
measured in the reference frame of the Earth. In Eq. [8j 
and pi± represent the momentum component parallel 
respectively perpendicular to the direction of motion of 
the center of mass. 

One more ingredient is needed to carefully analyze the 
proposed black hole signature, which is the angle be- 
tween the two showers in the Earth reference frame. One 
can start from (f> a and <f>b, which are the angles between 
the two emitted particles in the center of mass reference 
frame measured from the direction of motion of the cen- 
ter of mass (4> a + <j>b = it since the particles are moving 
back-to-back) and perform a simple Lorentz transforma- 
tion. The resulting angles in the laboratory reference 
frame are 



tan Qi 



IBHPBH^ + IBH COS! 



(9) 



These are the angles between the secondary showers and 
the direction of motion of the center of mass and the 
angle between the two showers is their sum. 

This distinctive black hole decay signature can be ob- 
served if the secondary showers are separated by an an- 
gle large enough for the experiment to be able to resolve 
the event into two distinctive coincident showers. This 
signature was already studied in detail for showers pro- 
duced in the atmosphere in [29 ! . As it is also specified 
in the reference cited above, numerical simulations show 
that there is parameter space for large angular separa- 
tion between the two secondary showers. The plot in 
Fig. [T] shows the angle of separation between the two 
secondary showers measured in the experiment reference 
frame (8 a ) as a function of the corresponding angle in the 
center of mass reference frame (<f> a )- One notices that for 
a range of values of the angle <j> a covering about 0.4° 
(179.8° < (f> a < 180.2°), 9 a takes large enough values for 
the showers to be resolved into two separate ones. 

The image in Fig. [2] shows the simulation of a quan- 
tum black hole decay into a muon and a proton {BH — > 
H~ + p + ) in water. Strictly speaking the black hole is 
formed by the scattering between the neutrino and a par- 
ton from the the nuclei. However for the sake of the sim- 
ulation is it best to consider the nuclei instead of its con- 
stituents. Needless to say, for all charges to be conserved 
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FIG. 1. Angle 9 a as function of the angle (j> a for m a corre- 
sponding to: m n + — 139.5 Mev and m^- = 139.5 MeV blue 
solid line, m p + = 938.2 MeV and m M = 105.7 MeV long doted 
red line, m p + = 938.2 MeV and m T = 1.777 GeV small doted 
blue line, (decreasing from a larger possible angle for the low- 
est value of m a , to a lower possible one for the largest value 
of ma- 



particles, combined with the large Lorentz factor of the 
center of mass, leads to the several orders of magnitude 
difference between the two energies. An angle of 5° be- 
tween the two particles in the Earth reference frame was 
used for this simulation (experiments with large enough 
resolution capabilities can possibly detect even smaller 
angles). For the sake of illustration, the axes of the plot 
have different scales, the angle thus appears larger that it 
really is. For the angle considered here, after the muons 
travel a distance of about a kilometer the separation be- 
tween the single muon and the bunch of muons remaining 




it was assumed that the black hole formed as a result 
of the impact between a muon neutrino and a neutron 
{y^i + n — > BH). It needs to be pointed out that this 
is just one of the possible decay channels for quantum 
black holes and it was chosen because it is suitable for 
the detection of the back-to-back decay signature (we do 
not mean to imply that this is the only decay channel for 
which the signature can be discovered). The numerical 
simulations were performed using CORSIKA-6600-WI- 
0.9 [43], (COsmic Ray Simulations for KAscade, a soft- 
ware package developed to perform detailed simulation 
of extensive air showers initiated by high energy cosmic 
ray particles - a version of which was also extended to 
simulate particle showers which develop in ice or water) 
@1I15]. More specifically the QGSJET 01C model for 
hadronic interactions at ultra- high energies was used |46j . 

In the Earth reference frame the energies of the two 
black hole decay products are approximately 10 17 eV for 
the proton and 3 x 10 11 eV for the muon. The reason for 
the large difference between the two energies comes from 
the fact that the angle in the laboratory reference frame 
can only take large values when the angles in the center of 
mass frame are close to 0° respectively 180° (see Fig. [IJ. 
The resulting large but opposite values of pii for the two 



FIG. 2. CORSIKA simulation of a quantum black hole decay 
into a proton and a muon in water. The black hole is produced 
by a 10 17 eV neutrino which collides with a neutron has a 
mass of approximately 10 eV and decays immediately. In 
the Earth reference frame the muon has an energy of 3 x 10 
eV and the proton has an energy of approximately 10 17 eV. 
The angle between the two particles in the Earth reference 
frame is 5°. 



from the proton shower is about 100 meters. One can see 
that the muon track (represented in green) is very well 
separated from the bunch of muons (black muon tracks) 
resulting from the shower produced by the proton, there- 
fore the muons can be identified as coming from separate 
showers but having the same origin. When analyzing the 
numerical simulations one also notices that most of the 
proton shower dies out very fast with the only particles 
propagating beyond a few tens of meters being a bunch 
of energetic muons. The muons can be identified by the 
Cherenkov light which is emitted by the particles trav- 
eling through water or ice. The situation is very similar 
for events taking place in ice and this is obvious if when 
comparing the plots in Figs. [3] and [I] 



Figs. [3] and [4] show side by side the number of parti- 
cles in the showers produced in water (respectively ice) 
by 10 17 eV protons (plots on the left side) and 3 x 10 11 
eV muons (plots on the right side). The figures on the 
left show that the only remaining components of the 
proton showers beyond some 30 meters are muons and 
Cherenkov light; while the hadronic and electromagnetic 



parts die out. The muons travel to distances of the order 
of kilometers while continuously generating Cherenkov 
radiation which can be measured by the underwater de- 
tectors, respectively by the detectors placed in ice. The 
figures on the right show the Cherenkov light produced 
by the single muons, which are in each case one of the 
two products resulting from the back-to-back decay of 
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FIG. 3. Plots representing the number of particles in the showers as a function of the distance traveled through water. Each 
plot contains ten overlapping numerical simulations. The left panel shows the particles and Cherenkov light produced in a 10 17 
eV proton shower, while the right panel shows the Cherenkov light produced by the 3 x 10 11 eV muon. 
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FIG. 4. Plots representing the number of particles in the showers as a function of the distance traveled through ice. Each plot 
contains ten overlapping numerical simulations. The left panel shows the particles and Cherenkov light produced in a 10 17 eV 
proton shower, while the right panel shows the Cherenkov light produced by the 3 x 10 11 eV muon. 



the quantum black holes. A 3 x 10 11 eV muon travels 
a distance of the order of one kilometer in water (simi- 
lar distances are traveled by muons moving through ice) 
and its track can be detected by the Cherenkov light 
which it produces. The muons remaining from the 10 
eV proton travel distances up to about two kilometers 
both in ice or water (the plots in Fig. [3] and [4] only show 
the shower development for the first five hundred meters 
in order for the electromagnetic component to be visi- 
ble). Beyond around on hundred meters, the number of 
Cherenkov photons emitted by the 3x 10 11 eV muons is of 
the same order of magnitude as the number of Cherenkov 
photons emitted by the muons generated in the 10 17 eV 
proton showers. 

As mentioned earlier, the original muon and the proton 
which generates the several other muons have very par- 
ticular energies for a certain initial neutrino energy. This 
is a direct consequence of the small range of angles in the 



center of mass reference frame for which the angle in the 
laboratory reference frame can have large enough values 
for the muon tracks to be spatially separated. When an- 
alyzing an event of the type described above one can also 
use the reconstructed energies to pinpoint back-to-back 
quantum black hole decay events. The price to pay is 
that only about 0.11% of the total number of quantum 
black hole decay events can be measured in this way. 

Measuring muon tracks with the characteristics de- 
scribed above in the data collected by a neutrino ob- 
servatory such as IceCube or the future KM3NeT [47] 
would signal a back-to-back black hole decay event. 

Factoring in the percentage of the total number of 
events which can be observed ( 0.11%) and the branch- 
ing ratios one expects the following number of events per 
year: 1, 21, 81, 179, 311, 440, 623, 806 for the first model 
of neutrino flux and 0, 2, 8, 17, 29, 40, 55, 73 for the 
second one (n = 0, 1, 2, 3, 4, 5, 6, 7). 
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IV. CONCLUSIONS 

We have described a novel approach to probe the pos- 
sibility for the Plank mass if it is in the few TeV range 
via back-to-back decays of Planck scale quantum black 
holes using the data collected by neutrino observatories. 

This signature consists in detecting the Cherenkov ra- 
diation of two simultaneous muon tracks (to be more 
precise one of the tracks is made up of several muons 
while the other would be from a single original muon) 
oriented at an angle, but pointing to a common origin. 
Another particularity which makes this signature unique 
is the energy in the laboratory reference frame of the two 
products of the back-to-back decay of the black hole. In 
the case studied here a 10 17 eV muon neutrino produces 
a 10 13 eV black hole by colliding with a neutron. The 
quantum black hole can be discovered only when the two 
particles resulting from its back-to-back decay are almost 
aligned (in the center of mass reference frame) with the 
direction of motion of the center of mass. In this case the 



energies of the two decay products (muon and proton) in 
the Earth reference frame are approximately 10 17 eV and 
3 x 10 11 eV; and the resulting muons produce two distin- 
guishable Cherenkov light tracks at an angle of several 
degrees. 

The available parameter space for this signature to be 
discovered is very small (this, of course, depends on the 
high energy neutrino flux) , but it is a very unique signa- 
ture and therefore worth being considered. 

This black hole decay signature allows the neutrino 
observatories to join cosmic rays experiments and the 
LHC in the efforts to search for TeV scale micro black 
holes. 
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